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Abstract 

Asymptotic behaviour of eigenvalues and eigenfunctions of a stiff problem is de- 
scribed in the case of the fourth-order ordinary differential operator. Considering the 
stiffness coefficient that depends on a small parameter e and vanishes as e — > on a 
subinterval, we prove the existence of low and high frequency resonance vibrations. 
The low frequency vibrations admit the power series expansions on e but this method 
is not applicable to the description of high frequency vibrations. However, the non- 
classical asymptotics on e of the high frequency vibrations were constructed using the 
WKB method. 

MSG: Primary 34E20; Secondary 74K10 

Introduction and main results. 

Stiff vibrating systems belong to a class of systems with singularly perturbed potential 
energy. Stiff problems are known in particular as boundary value problems for differential 
equations with very contrasting values of coefficients in different sub-domains. They relate 
to modelling vibrations of elastic systems consisting of two (or more) materials with one of 
them being very stiff with respect to the other. For the first time the stiff problems were 
investigated by J. L. Lions [1]. 

However, it is also of interest to describe the asymptotic behaviour of spectral properties 
for the stiff problems. These system has two types of eigenvibrations, namely low frequency 
vibrations and high frequency ones. From a physical viewpoint we can postulate that two 
kinds of eigenvibrations can appear: one for the stiffer structure and the other for the 
softer structure. Different aspects of the spectral stiff problems are considered in [2]-[10] 
with the best general reference being [8]. The asymptotic behaviour of the low frequency 
vibrations has been widely studied with different techniques [2], [3], [7] and [9]. 

In this paper, following [10] we consider the phenomenon of high frequency vibrations. 
The leading terms of high frequency vibrations for different problems were constructed 
in [7-9]. Information on the behaviour of high frequency vibrations was also provided in 
[10]. Studying the stiff problem for the forth-order differential operator, we construct the 
complete asymptotic expansions of high frequency vibrations using the WKB technique. 
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Why do high frequency vibrations appear? On the one hand, the spectrum {A|}^j^ 
of the stiff problem is asymptotically dense in [0, oo) as e — (see [10]). On the other 
hand, for fixed i the asymptotic expansions of an eigenvalue Xf and the corresponding 
eigenfunction uf are nonuniform with respect to i. Hence the sequences u?^^^ with i{e) — > cxd 
can support stable forms of vibration as e ^ 0. In addition, the corresponding sequencies 
of eigenvalues A| converge to some positive limit points (see Fig. 1, 2). It is also shown 
that the approximation of the limit form of vibrations by sequence {ttf } has the discrete 
character. Therefore we construct and justify the asymptotics for a family of discrete sets 
of a small parameter. 
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Fig. 1 Fig. 2 

Moreover, such approximation is ambiguously determined. Then we introduce a defor- 
mation parameter 6 in the asymptotics. Consequently, for each 6 we construct expansions 
of u^, though they are asymptotically equivalent when e — > 0. 



1 Problem statement 

Let an interval (a, 6) C M contains the origin. Let us consider the eigenvalue problem 

Ue{a) = <(a) = 0, u,{b) = <(6) = 0, (2) 

Ue{-Q) = Ue{+Q), <(-0) = <(+0), (3) 

(fco(x)<)(-0) = e4(fci(x)<)(+0), (A:o(x)<)'(-0) = e\k^{x)u'i)' {+^) . (4) 

Here e is a small positive parameter and a function 

kaix), a < X <Q 
e'^ki{x), < X <b 



ke{x) 



is smooth and strictly positive in (a, 0) U (0, b). We study the asymptotic behaviour of the 
eigenvalues and the eigenfunctions of (1)~(4) as e ^ 0. 
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The problem models eigenvibrations of a non-homogeneous rod. The rod consists of 
two components with the same density of mass and strongly different elastic properties. 
The fourth power of e in the definition of fcg is suitable for the next consideration (see 
Section 4). 

Let us introduce the Sobolev space H^la, b) as the closure of set C^{a, h) with respect 
to the norm 

,b \ 1/2 

\u\\ = I / dx 



and bilinear forms in HQ{a, h) 

i-O rb 

ao{(p,'ilj) = / koip"'il^" dx, ai((/p, V')= / knp"'il;" dx, ae = ao + e^ai. 

Ja Jo 

For each e > we note by || • ||e the norm in HQ{a, b) associated with the form «£(•, •). It 
is known as an energetic norm associated with the elastic-like energy. 

Let us consider a variational formulation of (l)-(4) that is to find Ag and € Hq^u, b) 
such that 

as{Ue, if) - Xs{Ue, ^)L2{a,b) =0, if e ifg («' (5) 

Problem (5) is a standard eigenvalue problem with a real discrete spectrum. For each fixed 
e let us consider an eigenvalue sequence 

< Af < A| < • • • < Af < . . . , where Af ^ oo as i ^ oo. 

Note that each eigenvalue is simple. Let the corresponding eigenfunctions {u|}^;^ form an 
orthonormal basis in L2{a, b). 

2 Asymptotic expansions of low frequency vibrations 

We investigate the asymptotic behaviour of eigenvalues Af and eigenfunctions uf for a fixed 
number i. 

Lemma 2.1. Each eigenvalue A| is a continuous function with respect to the parameter e, 
e £ (0, 1). Moreover, 

Af ^ Cie\ 

where constant Ci is independent of e. 

Proof. The continuity of eigenvalues follows from the variational principle 

Af = sup mf —2 , (6) 

where P is a (i— l)-dimensional subspace of i?Q (a, b) and P"*" is the orthonormal complement 
of P. 
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Suppose that the supremum m (6) is achieved on a subspace Pj. Since Pi is finite 
dimensional, we choose a function /j € \ {0} that vanishes in (a, 0). Then 

^ ao(/, /)+£%(/,/) 

^ aoifi, fi) + £^ai{fi, fi) _ 4 aiifiJi) _ ^ 4 

||f.||2 - ^ ||#.||2 - ' 

since ai(/i,/i) / 0. □ 

We postulate the expansions of the eigenvalue = A| and the eigenfunction = «| 
for given i € N: 

A, ~ e4(Ao + e^Ai + ...), 

e^{uQ{x) + e^ui{x) + . . xe{a,Q), (7) 
■uo(2;) + e''t;i(x) + . . . , xG(0, 



Ue{x) 



(8) 



Substituting (7) into (l)-(4) we deduce that Aq is an eigenvalue and vq is an eigenfunction 
of the problem 

{kiVQ)"{x) - \qvo{x) = 0, X e (0, h), 

7;o(0) =^;o(0) = 0, vo{b) = v'^ib) = 

Note that all eigenvalues of the problem are simple. 

Then the function uq is a solution of the boundary value problem 

(A:oO"(x) = 0, xe(a,0), 
no(a) = u'^{a) = 0, A;o(OX(0) = A:i(0)<(0), (A;oO'(0) = (A:i^;^;)'(0). 

Let us define next terms of expansions (7). The function vi satisfies 

{kiv'l)"{x) - Ao'Ui(a:;) = \iVq{x), x e (0, 6), 
vi{id) = t/o(0), v'M = ^^o(0), vi{h) = v'^{h) = 0. 

Since Aq is an eigenvalue of (8), such problem has a solution under the condition 

^1 = 11^0 11^2(0,6) ((^i^o)'^^o - kiv'^u'o) 
Further, the function ui is a solution of 

{kQu'{)"{x) = \qUq{x), X G (a,0), 
ni(a) = <(a) = 0, A:o(0)<(0) = A:i(0)<(0), (A;o<)'(0) = (A:i<)'(0). 

The general terms of (7) can be found in the same way. The asymptotic expansions 
are justified by classical methods [14, 2]. 

The vibrations described above are corresponding to the low levels of potential energy, 
since a^iue, Ue) = o(e^) as e ^ 0. Naturally we have vibrations of the soft part of a system 



x=0 
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only. The soft part is clamped at x = by the immovable stiffer structure. Note that the 
leading terms of expansions (7) are determined by the soft part. Since every eigenfunction 
uf converges to an eigenfunction vq of (8) extended by zero to (a, b), the set {uf}^^ is not 
a basis in L2{a,b) for e = 0. In other words, the "low frequency" region does not provide 
a good insight on the vibration problem over all (a, b). Therefore we consider here another 
kind of vibrations, namely, vibrations with "finite non- vanishing energy". 

3 High frequency vibrations 

On the one hand the asymptotics of and u| are nonuniform with respect to i, on the 
other hand the spectrum of (l)-(4) is asymptotically dense in the positive spectral semi- 
axis (see Lemma 2). On account of the above remarks, we can find stabile vibrations ti^^^^ 
associated with certain sequences of eigenvalues A|^^^ with i{e) oo. 

Let us denote by E the set of all pairs (e, A(e)), where A(e) is an eigenvalue of (1)"(4) 
for some e G (0, 1). 

Lemma 3.1. The closure of the set E includes semi-axis {(e. A) : e = 0, A ^ 0}. That is 
to say each positive A can be approximated by a sequence of eigenvalues A(e). 

Proof. Suppose that for a certain A > there exists a neighborhood U of the point (0, A) 
such that U D E = 0. Let us choose e* > small enough for (e*. A) € U. Recall that Aj(e) 
is a continuous function with respect to e and Aj(e) — > as e ^ (see Fig. 1). Since the 
neighborhood U by definition does not contain points (e* , Xi{e*)) then Aj(e*) ^ A for all 
i G N that is impossible. □ 

Let us consider a sequence of pairs (ej,Aj) G E that converges in to (0,u;^) with 
Lo > as i ^ oo. Let Ui be an eigenfunction of (l)-(4), being normalized in ^2(0, 6) and 
associated with Aj. 

Definition 3.1. We say that the eigenfunction sequence {ui}"?^^ supports high frequency 
vibrations asi ^ 00, if the sequence of restrictions g) ^ nonzero limit v in H^{a, 0). 
Then uj is referred to as the limit frequency and v as the limit form of these vibrations. 

Remark 3.1. It will be shown that the functions Ui have strongly oscillatory character in 
(0,6). Moreover, the vibrations ui have the "energy" close to uj'^ as i ^ 00. 

Theorem 3.1. // a sequence {ui}^^ supports high frequency vibrations with a limit fre- 
quency UJ and a limit form v, then uj'^ is an eigenvalue and v is an eigenfunction of the 
problem 

{kQ{x)v") -Lv^v = X e (a,0), 
v{a) = v'{a) = 0, v"{0) = v"'{0) =0. ^ ^ 

and restrictions tij|(o.6) converge to zero in the weak topology 0/^2(0,6). 
Proof. From (5) we have ||Mi||e = Ai||uj|| 2,2(0,6) = Aj. Then 

e'^ai{ui,ip) ^ e'^ay^ {ui,Ui)ay'^ {if, if) ^ \\ui\\e^ay'^{ip,ip) ^ Aia^^^((/3, v?) (10) 
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for any G HQ{a, b). Let us consider identity (5) for test functions if € C^(0, b) 
Using (10) we obtain 

It follows immediately that Uj ^ in L2(0, b) weakly. 

Throughout the proof, Ha denotes the space {v G H^{a,0) : v{a) = v'{a) = 0}. Passing 
to the limit as — > in (5) we obtain 

ao{v,ip) - X{v,ip)L^(^a,o) =0, ipeHQ{a,b). (11) 

Since the restriction (/'|(o,a) is an element of Tia, the identity (11) corresponds to the eigen- 
value problem (9). This finishes the proof, because v is a nonzero function. □ 

Consequently, the limit frequencies uj are generated by means of stiffer part of a vibrat- 
ing system. 

4 Asymptotic expansions of high frequency vibrations. The 
leading terms 

By We we denote an eigenfrequency of vibrating system (l)-(4). We search for lOs with 
the asymptotic expansion 



J^e'^u;,., LOoy^O. (12) 



oo 
A;=0 

Let us consider a one-parameter family of the following vector-functions 



N{>c,x) = (cosx5(x), sinx5(x), e''^^^^'), e^(^(^')-^(^))) , x G (0,&), x € M, 

where ^ 

S{x) =uJo [ k~^'^{t)dt. 





Let < •, ■ > be an inner product in M^. We postulate the expansions of in the form 

/ oo 



Ue{x) ~ < 



Y.e'^Vkix), xG(a,0), 

fc=0 

oo 

Y^e^ <h{x),N{e-\x)>, x€(0,6), 
k fe=o 



(13) 



with A : (0,&) ^Ml 

Remark 4.1. On the interval (0,6) we have the equation 
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with a small parameter at the highest order derivative. Thus the construction of expansions 
(13) on the interval (0,6) was motivated by [13] with the method of WKB-approximations, 
which are also known as short-wave approximations. Hence, we can find a solution in 
the form 



oo 



e ^ flfcW 



fc=0 

where function S{x) has to satisfy the eikonal equation ki{x)S' (x)'^ — uJq = 0. Since there 
exist 4 different solutions of the eikonal equation with respect to S' , we have introduced the 
function N{k, •) combining in it's components a fundamental set of solutions. 

Substituting (13) for in conditions (4), we obtain 

feo(O) v'^iO) + ...^eHi{0) 5' (0)2 < /o(0), T^N {e~\o) >+..., (14) 
{kov'^y{0) + ...r^eki{0)S'{Of <fo{0),T^N{e~\0)>+... . (15) 

Above we take into account the equahty 

■^N (e-\x) = e-^S'{x) TN (e-\x) 
dx 

with an orthonormal matrix 

^=(^' rj' ^i = (? "o^)' ^2 = (V \ 

By (14) and (15), we obtain ^^(O) = v'^' {^) = 0. Substituting expansions (12), (13) 
for iOf, and in equation (1) and boundary value conditions (2), we have the eigenvalue 
problem 

d'^ ( d'^vo\ 4 / N 

ko{x)-j-Y ) -ujqVo{x) =0, xe{a,0), 



dx-' V dx' J ^° ' " " ^ (16) 

vo{a) = v'oia) = 0, v'^{0) = v"'{0) = 0. 

Let ujQ and vq be an eigenfrequency and an eigenfunction of the problem (compare with 
Th.l). Note that every eigenvalue loq is simple. Suppose that the function vq satisfies 
condition \\vo\\L2{a,o) = 1- 

In order to find the leading term /o of expansion (13), we substitute the second of series 
(13) in equation (1). In particular, we obtain 

(fci5'^-a;4)/o = 0, (17) 
AkiS'^Tf'^ + {2k[S'^T + QkiS'^S"T - 4a;|]a;i/)/o = -{kiS'^ - uj^fi. (18) 

Since the function 5 satisfies the eikonal equation kiS'^ — = then (17) holds. 
Hence (18) is a homogeneous system of linear differential equations with respect to the 
function /q. 

Taking into account conditions (3) and boundary value conditions (2) at x = 6, we can 
write 

f'^ = A{x)fo, xG(0,6), 
(/o(0), iV(e-\ 0)) = 7;o(0), (/o(0), TN{e-'\ 0)) = 0, (19) 
(/o(6), N{e-\h)) = 0, (/o(6), TN{e~\b)) = 0. 
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with a matrix 



A 



\ 



^1 






^1 

k' 











■ — 

1 





8fci ^ ^ / 

Note that the matrix A depends on a parameter lvi. We shah define wi below. 

On the one hand problem (19) depends on e by means of the boundary value conditions, 
on the other hand that one is an ill-posed problem. To resolve both these problems, we 
consider a discrete set of small parameter e. We shall choose below a specific sequence 
— > of a small parameter and asymptotic expansions will have a discrete character with 
respect to e. That agrees with a discrete phenomenon of high frequency vibrations. 

Lemma 4.1. Let w : (0,6) — > M"^ 6e a smooth vector-function and a be a vector in M^. 
There exists a infinitely small sequence {ep}^i such that the problem 

y'{e,x) = A{x)y{e,x) +w{x), xe(0,6), (20) 
{y{e, 0), N{e~\0)) = ai, {y{e, 0), TN{e-\Q)) = ^2, (21) 
{y{e,b),N{e-\b)) = a^i, {y{e,b),TN{e-\b)) = a^, (22) 

has a unique solution y{ep, •) for p € N. The family of solutions {y{£p, Olp^i holds the 
inequality 

_M_ 

Wvi^p, •) - y*llci «S C'e 

for a smooth function : [0,b] M'^. The constants C and M are independent of 
parameter e. 



Proof. The fundamental matrix of (20) is 



/ cos^5(x) 



$(x) = k-^l^{x) 



sm'^Six) 



COS "^Six) 



\ 
















"0 ^ ' 










5l(5(x)-5(fe)) 



Therefore we have a representation of the general solution 

y[x) = ^{x)[l3^h(x)), 

where /3 is a constant vector and h{x) = ^~^{t)w{t) dt. Suppose a vector-function 
y{e,x) = <I>(x)(/?£ -|- h{x)) is a solution of (20)-(22). Seting 7^ = +ujiuJq^, it is easy to 
check that 

^\x)N{e-\x) = k^^^^{x)N{js,x), (23) 
where is a transposed matrix. Then 

< yi£,x),N{e~\x) >= k^^^\x) < (3e + h{x),N{-ie,x) >, 
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and we can write (21), (22) in the form 

iV(7.,0)) = 
(/3„riV(7.,0)) = 
N{j„b)) = 

{(3e,TN{-f„b)) = 



modi, 
moo-2, 

micjs- < h{b), N{je,b))), 
miai-{h{b), TiV(7„6))), 



(24) 



where tuq = A;J^^(0) and mi = ky^{b). Note that the matrix <I>* commutes with T and, 
moreover, h{0) = 0. 

Hence, the vector /3e is a solution of a hnear algebraic system with a matrix 



G(7e 



/ 1 

1 

cos7£S'(6) sin7eS'(6) 

V -sin7£5(6) cos7e5'(6) 



-7.5(6) 



e-7.S(6) \ 
e-7.5(6) 

1 
1 



e-^^^^'') cos 



leS{b)) 



Since S{b) ^ 0, the determinant 

detG(7s) = -2 cos 7^5(6) + 2e-T=^('')(2 

does not vanish for all e > 0. 

Let us fix 5 € [0,27r) and choose a sequence Ep from the set of conditions ^epS{b) 
5 + 2Trp for p = 1,2, . . . . Namely, let 

ujoS{b) 



oJo{6 + 2Trp)-uJiS{b) ^^^^ 

for all p ^ pq. Here po is the smallest natural number such that the denominator of (25) 
is positive. Without restriction of generality we set po = 1. We denote 7p = jsp- Since 
7p — cxD as Ep — s- 0, the matrix G{^p) is an exponentially small perturbation of a matrix 



Gn 



1 
1 

cos 6 sin 6 

- sin 5 cos 5 



1 0\ 

-1 

1 

1 / 



Since N{jp, b) = (cos (5, sin 6, e >'^(''), then the right-hand side of (24) is an exponen- 
tially small perturbation of vector 

g = (mocTi, mofT2, mias - {h{b), Ns), mia^ - {h{b), TNs)), 

where Ng differs from vector N{^p, b) by the third component only. 

Let us suppose that 6 is not equal to tt/2 and 37r/2. Then the matrix Gq is non- 
degenerate. From the theory of finite-dimensional perturbations we obtain 

||/3e,-/3.||M4^Ce->^(^), 

where is a solution of Go/3 = g. Let y*(x) = <I>(x)(/?* -|- h{x)), then 

Wvii^p,-) -y*\\c' ^ ll^llci|l/5ep - /3*IIr4, (26) 
<l>(x)|| -|- ||$'(x)||). Note that 7^ ^ Coe""*^ with a positive constant 



where II^Uci = max 
xe(o,b) 

Co, which proves the Lemma 



□ 



9 



Remark 4.2. From now on, we say that y^, is a solution of (20)~(22) with neglect of 
exponentially small terms. However the choice of a sequence £p is non-unique and depends 
on 6 and uji . We shall define u)i at the next step hut, on the other hand, we shall keep the 
dependence on 6, which will be a deformation parameter. The approximation of the limit 
function vq by eigenf unctions Un{e) ambiguously determined. Hence we can not define 
5 uniquely. Consequently, for each 5 we construct the expansions of , though they are 
asymptotically equivalent as e ^ 0. 

Let us return to the study of problem (19), that is a sub-case of (20)-(22) with the 
right-hand side vu = and a = (t;o(0), 0, 0, 0). Since (19) is a homogeneous system, we 

obtain /o = $/3oi where Pq = hk]^^ (0)vo{0) ( 1 — tan J, 1 -|- tan (5, 1 -|- tan 5, ) is a 

\ cos J 

solution of the corresponding system with matrix Go- 

In order to calculate the first-order correction uji, we consider the problem for vi 

^ (^^o(a;)^^^ - ujqVi = Aujiuj^vo, x G (a, 0), 

k^^^(O) ^^^^ 
v,{a)=v[{a) = 0, v';{0) = 0, <(0) = -4-^u;o%(0)(l + tan 5). 

Since uJq is a simple eigenvalue of (16), there exists a solution of (27) under the condition 
ui = —jky^{0)v'^{0){l +tan6). Let us choose a solution vi such that (t'l, fo)L2(a,o) = 0- 
Now by (25), it follows that 

eA^) = , (28) 

4cJo(5 + 27rp) + A:|/^(0)t;g(0)5(6)(l +tan5) 

where 5 € [0, 27r) \ {vr/2, 37r/2}. Let us denote by £s the sequence (28) for fixed 5. We shall 
construct the asymptotic expansions only for e E £s. Hence, we have just found loq, vq, /q, 
tgcorrections uji, vi and the sequence £s- Recall that /o, wi and vi depend on 6. 



5 Complete asymptotics of high frequency vibrations 

Let us find the general terms fk, iOk+i and Vk+i of expansions (12), (13) for k ^ 1. In the 
same way as for the vector /o, we obtain the boundary value problem for /fc. With neglect 
of the exponentially small terms in the boundary conditions, we can write the problem in 
the form 

fl =A{x)fk + Wk, xe{0,b), 
(/fc(0), N{e-\0)) =VkiO), {fk{h), N{e-\b)) = 0, 
(A(0), TN{e-\0)) =u^\kl^\o)v',_,{0)-kl/'{0){'^-HO)fU{0)^No)), 
{f,{b), TN{e~\b)) = -u^'kl/\b){^-Hb)f,_,{b),Ns), 

where Nq = (1, 0, 1, 0) and the A'^^ is defined in the proof of Lemma 3. 

For each k the right-hand side Wk of (29) is a smooth function, namely, 

1 ''"^ d 
Wk = ^^(E >^k^^T^fm - S'T^Pk - -^ik[{S'fTfk^i+ 

^ m=0 



3kiS'S"Tfk-i + 3fci(50'T/^_i + T^Pk-i)) 



(29) 
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where Pi = hS'T'^f'/^^ + £{2kiS'T^f^_^ + kiS"T^fi^i + kif/_^). The proof is by induction 
on k. According to Lemma 3 there exists a solution of (29) for e G £s. 
Now we can find Wfc+i and Vk+v 



^(k,{x)'^\ 



k+l 



^j,2 j '^U'^K+i — ^ ^ ^mVk—m+l-i X € (o, 0), 

^ 771=1 

t;fc+i(a) = 0, v'f^^^{a) = 0, 
(«+i)(0) = {<^-H^'okfT'fk^, + k-'/'Q,^2),No) 

{kov'^iYm = {^~Hu;lky'Tf, + k-'/'Rk-i),No) 



(30) 



x=0 
x=o' 



where 



Qk = h — iS'T^fk + fU) + S'T^f'k, 
ax 



d 



S'T'Qk-i + j^{ulk\'^T^fk + Qfc_i) 



and Am = X] ^i^j^i^s with i + j + / + s = m. Since is an eigenvalue of (16), boundary 
value problem (30) has no solution for an arbitrary right-hand side. We can write the 
existence condition in the form 

v'Q^-\iolk\'''Th + A:-'/'i?fc_i), A^o > 

Let us choose a solution Vk+i of (30) such that (wfc+i, 't^o)L2(a,o) = 0- 
Hence, the algorithm scheme of asymptotics (12), (13) is 

loq ^ vq ^ fa ^ uji ^ £s ^ vi ^ > fk^ Wfc+i Vk+i - 

Note that all terms, except for cjq and vq, depend on parameter 5. 



x=0 



6 Justification of asymptotics 

The nonstandard object of investigation, namely, the sequences of eigenfunctions u^i(^^y 
changes the classical scheme of justification. Note that in Sections 4 and 5 we have con- 
structed series (12), (13), however we have not defined the object that is approximated by 
ones yet. 

Only by using the formal series (12) we shall define a sequence u^^^-^ that supports the 
high frequency vibrations with the limit frequency ujq and the limit form vq. 

Let us fix n € N and introduce a sequence of real numbers {Ae^^jggf^ and sequence of 
functions in H'^{a,b). Namely, for each e E f^^ we set 

A(") = Ao + eAi -h . . . e"A„,, A^ = ^ LOiUjjUJkUJi, (31) 

i+j+k+l=s 

vo{x) + evi{x) -\ \-e'^Vn{x), x e {a,0), 

4"^(^) = <| f2e^<f,{x),N{e-\x)>, x e {0,b), ^^^^ 

i=0 
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where the numbers ujk, the functions v^, the vectors and the set £s are defined in Sections 
4 and 5. 

Lemma 6.1. There exists a sequence {X'^}ee£g of eigenvalues for (l)-(4) such that 

|AW-A^| ^C„e"+\ £G£s- (33) 

Proof. Denoting by He the space H'^{a, h) with the inner product as{-, •), let us introduce 
an operator such that 

Note that is self-adjoint and compact for all e > 0. Then we can write problem (l)-(4) 
in the form 

AeUe - {X')-^Us = 0. 

Substituting sequences (31), (32) in problem (l)-(4), we obtain 



n(")(a) = 0, -;^nW(a) = 0, 4")(6) = 5l^)(e), ^u^J^Hb) = g^^He), 



d_ 

dx 



— ( ke{x)^ut^ ] - A(")4") = F„(e,x), X G (a, 6), 



dx"^ 



,in) 



-u. 



(n) 



d_ 

dx 



dx 

dx^ ^ dx"^ ^ ^ 



(34) 



where [/]o is a jump of a function / at x = 0. The right-hand sides of problem (34) satisfy 
the inequalities 



\Fn{e,x)\\co^a,b) ^ |#(e)| ^ C^e"^, i = 1,2,3, 

\hn{e)\ ^ C„e", \z^\e)\ ^ C„e"+\ i = 1,2. 



(35) 



The function ul does not belong to the space Tie, because it has a point of discontinuity 



at X = 0. Taking into account (35) we can choose a function (^^"^ such that u^'^' + ipe"'' S Tl, 



and 



max 

x6(a,ti) 



4") 



+ 



+ 



^2 



(36) 



Let Ve^"^ = >Ce{u^J^'^ +Lpi'^), where is a normalizing constant such that ||V^''"'* H-h^ = 1- It 
is easy to check that ^ xq > 0. It follows from (35), (36) that 



[a, - {xt^r^i) 



(37) 



where is a constant independent of e. Hence, according to the Vishik-Lusternik lemma 
[14] there exists the eigenvalue (A^)~^ of the operator A^ such that 



1 1 

A^ 



A^) 
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Applying this inequality for the value n + 1 instead of n, we obtain 

□ 

Lemma 6.2. There exists a positive number d such that for each e ^ £s the interval Is = 
(A^"^ — de,Ai"^ + ds) contains exactly one eigenvalue X^^^^ of problem (1)~(4)- Moreover, 
the eigenvalue number k{e) satisfies the inequality 

ao£-^ ^ k{£) s: aie~^, (38) 

where oq, oi are constants independent of e. 

Proof. According to Section 2, we have 

A^ = eVm(l + am(e)), £ ^ 0, 

where fim is an eigenvalue of (8) and = o(l) as e ^ 0. Moreover, the following 

asymptotics hold 

fim = "T-^lco + T{m)) with r(m) = o(l) as m ^ oo. 

Then 

A^ = (em)^(co + r(m))(l + a„,(e)). (39) 
By Lemma 4, for each e ^ £s there exists a number /c = k{£) such that 

Ao - h^e ^ A|(^) ^ Ao + hoe 

with 6o > 0. We conclude from (39) that 

Ao - feoe ^ {ekf{cG + r(A:))(l + afc(e)) ^ Aq + h^e, 

and finally (38) holds. According to (39) we have 

Xl^^^ = coe^k{ef + o{e^), Ss^e^Q. 

Then the distance between the two neighboring eigenvalues that are the closest to the point 
Ao admits the estimate 

|A|(,)+i - A|(,)| = ^c,e^k{ef + oie'), £s 3 e ^ 0. (40) 

Taking into account the asymptotic behaviour of k{e), we obtain the existence of the 
interval Ig of the length 2d£ such that exactly one eigenvalue A^ of problem (l)-(4) belongs 

(n) 

to it. The case of the point A^ being a midpoint of is impossible since this contradicts 
(33) and (40). □ 

Now we can introduce the object that is approximated by formal series (13) constructed 
in Sections 4 and 5. Let {A|^^^}£g£:^ be the sequence which is defined in Lemma 5, where 

A|^^-j is the nearest eigenvalue to the value Ae"\ Let {ul^^^j^g^-^ be the sequence of the 
corresponding eigenfunctions, Hw^ = 1- 



13 



Theorem 6.1. As £s ^ ^ ^ the sequence {u^f,(^^-j}ee£s supports the high frequency vibra- 
tions with the limit frequency loq and the limit form vq . Moreover, 

eeSs, (41) 

for n = 0, 1, . . . . In particular, the following inequalities hold 



En 
Vkix)e' 

- ^eY^l {fk{x),N{e-\x))e'' 



(f 

where = HV^^^^Hr^ is a normalizing constant and ^ xq > 0. 



Ci(0,fc) 



n-1 



(42) 



Proof. By Lemma 7, we choose do > such that do^-neighborhood of the point (A^"^) ^ 
contains exactly one eigenvalue {^k{e)^~^ °^ operator A^. From inequality (37), we 
obtain (see [14]) 

||n,-V;(")||, ^i^„do'^""'' (43) 
where is a normalized eigenfunction of A^ associated with the eigenvalue (Al^^-^)"^. 
Hence, = zbit^j.^^. There is no loss of generality in assumption that (43) holds for the 
function u^^y Taking into account that the terms of series (13) are bounded in T^g-norm, 
from (43) for the number n + 2, we obtain (41). 

In addition, by (36) and the definition of norm, we obtain 

\\Uk{e) - ^eUl \\m{a,0) ^ We , |t%(e) " XeU\ '\\m{Q,b) ^ We 

Estimate (42) follows from the Sobolev embedding theorem. 

It is easy to check that the constant has a nonzero limit as e — > 0. Then sequence 
'L'"fc(e)J'ee£'i Converges to the eigenfunction kvq of problem (8) in (a, 0), namely, this se- 
quence supports the high frequency vibrations. □ 
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